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Abstract 

When  ocean  waves  in  deep  water  interact  with  a  current,  the  direction  of  propagation 
andcharacterist.es  of  the  waves  snch  as  height  and  length  are  affected.  SJlTZ£ 

the  Gulf  Stream  or  Antarctic  Circumpolar  Current.  Remote  sensing  techniques  such  as 
synthetic  aperture  radars  (SAR)  have  the  potential  to  detect  wave  systeTover  tid" 
wSp  ""•  COmbit6  a  m°dd  fOT  """  refrMti0n  »  "»  P-—  2  ILJnU 
tudv  the  t™    ^  ^C"0n  aDd  magnitude  °f  the  ^-vening  cmrents.   In  this 

e^leo  T?or  1 7  Tl  by  SAR  on  a  sateUite  pass  over  the  Gulf  st— 

urr^T  fi,  ^T^?  PIedjCted  *»  '  nUmeriCaI  modd  ™"«  «»*«<»•  <*  varying 
h^uiH  T  T:  T"eS  ',S  C°mPaIed  '°  SAR  geaerated  ™e  *«*»•  By  mSnl 
dl^Tf  P  ^  ir"145  fa  *he  bKt  ""^tion  of  wave  refraction  to  the  SAR 

data,  the  tomographic  problem  of  measuring  the  Gulf  Stream  current  is  solved 
i,    ~^Tl  Cm^a  be'ween  the  mod«l  ««»  SAR  data  i.  obtained  when  a  current 
mt"  bJ.  a   °P  hat  TdocitT  ?">««  «*>  »  direction  of  75'  and  a  current  speed  of 
2  m/..  The  direction  agrees  with  that  visually  observed  from  the  SAR  image.,  and  the 

ttfe  o'^^Tclt  "*  *£  '°  ^  COaS'  Guard  eStimates  for  the  Gulf  Stream  at  the 
time  of  the  SEASAT  pas..  The  current  profile.  used  did  not  take  into  account  a  po.sible 
widening  of  the  Gulf  Stream  at  the  position  of  the  sateUite  overpa...  There  i.  a  (treat 
deal  of  .catter  in  the  SAR  data,  both  before  and  in  the  Gulf  Stream,  „  it  iTafficuU  t 
correlate  every  point  with  .pecific  current  behavior,  but  the  increase  in  wave  length  and 
change  m  wave  angle  m  the  center  of  the  Gulf  Stream  seem  to  indicate  that  there  may 
be  a  nouWorm  feature  such  a.  the  formation  of  an  eddy  or  other  lateral  variability 
near  the  current  s  edge.  J 
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Chapter  1 
Introduction 


1.1      Motivation  and  Problem  Statement 

When  ocean  waves  in  deep  water  interact  with  a  current,  the  direction  of  propagation 
and  characteristics  of  the  waves  such  as  height  and  length  are  affected.  The  influence  of 
currents  on  wave  propagation  is  more  pronounced  in  deep  water,  because  as  the  waves 
shoal,  depth  refraction  dominates  the  changing  wave  characteristics.  Swell  in  the  open 
ocean  can  undergo  significant  refraction  as  it  passes  through  major  current  systems  like 
the  Gulf  Stream  or  Antarctic  Circumpolar  Current.  Remote  sensing  techniques  such  as 
synthetic  aperture  radars  (SAR)  have  the  potential  to  detect  wave  systems  over  a  wide 
geographical  area.  Combining  a  model  for  wave  refraction  in  the  presence  of  currents, 
with  SAR  measurements,  the  inverse  problem  of  using  the  measured  wave  data  can  be 
solved  to  determine  the  direction  and  magnitude  of  the  intervening  currents.  In  this 
study  the  behavior  of  swell  measured  by  SA_R  on  a  satellite  pass  over  the  Gulf  Stream 
is  examined.  The  refraction  predicted  by  a  numerical  model  under  conditions  of  varying 
current  profiles  and  velocities  is  compared  to  SAR  generated  wave  spectra.  Complex 
correlations  are  calculated  between  the  vectors  formed  by  SAR  measured  wave  length 
and  direction  and  the  wave  length  and  direction  vectors  predicted  by  the  numerical 
model.  By  chosing  the  model  inputs  which  result  in  the  best  correlation,  the  tomographic 
problem  of  measuring  the  Gulf  Stream  current  is  solved. 
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In  a  pioneering  work,  Unna  (1942)  showed  how  waves  are  shortened  or  lengthened 
by  opposing  or  following  currents,  respectively.  Johnson  (1947)  derived  a  form  of  Snell's 
law  which  describes  the  refraction  of  wave  crests  when  interacting  with  a  current  at  an 
oblique  angle  to  the  wave  propagation  direction.  Arthur  (1950)  introduced  the  concept  of 
wave  energy  transmission  along  rays  which  are  not  necessarily  perpendicular  to  the  wave 
crests.  Longuet-Higgins  and  Stewart  (1960)  demonstrated  the  importance  of  nonlinear 
interactions  between  waves  and  currents  and  introduced  the  concept  of  radiation  stress  to 
describe  the  variation  of  wave  energy  in  the  presence  of  a  moving  medium.  Kenyon  (1971) 
derived  an  analytical  solution  for  the  refraction  of  orthogonals  and  rays  in  the  presence 
of  a  linear  current  shear.  Abemethy  and  Gilbert  (1975)  developed  a  numerical  model  to 
predict  the  refraction  of  surface  gravity  waves  over  an  ocean  bottom  of  varying  depth. 
Treloar  (1985)  modified  this  model  to  include  the  linear  effects  of  current  interactions. 

A  depth  refraction  model  developed  by  Brampton  (1977)  is  modified  to  solve  the 
kinematic  problem  of  wave  refraction  in  currents.  This  current  refraction  model  is  val- 
idated by  comparing  it  to  analytical  solutions  derived  by  Johnson  and  Kenyon.  It  is 
then  applied  to  an  ocean  situation,  and  the  wave  kinematics  observed  by  a  synthetic 
aperture  radar  flown  on  SEASAT  during  a  pass  over  the  Gulf  Stream  are  compared  to 
the  kinematic  behavior  predicted  by  the  model. 

1.2      Outline  of  Present  Study 

In  section  2  the  theory  of  linear  wave  interactions  with  currents  is  reviewed.  The 
governing  equations  of  wave  refraction  in  the  presence  of  currents  are  shown,  and  the 
derivations  of  analytical  solutions  to  the  refraction  of  wave  fronts  passing  over  a  disconti- 
nuity, and  waves  passing  into  a  linear  shear  current  are  reviewed.  Then  the  assumptions 
and  structure  of  the  numerical  model  used  in  this  study  are  described,  and  a  number 
of  examples  are  shown  comparing  analytical  and  numerical  solutions  to  the  problem  of 
wave  refraction. 

In  section  3  the  principles  of  imaging  real  waves  with  synthetic  aperture  radar  are  re- 
viewed, and  the  characteristics  of  the  instrument  as  deployed  on  SEASAT  are  described. 


The  SEASAT  data  from  pass  1339  used  in  this  study  are  discussed. 

In  section  4  the  application  of  SAR  data  and  the  numerical  model  for  predicting 
wave  refraction  in  the  presence  of  a  current  is  applied  to  the  Gulf  Stream.  First  the 
environmental  and  meteorological  conditions  which  existed  at  the  time  of  pass  1339  are 
presented.  Next  the  analysis  of  SAR  images  to  produce  wave  energy  spectra  is  described. 
Using  these  spectra,  a  wave  system's  propagation  is  traced  across  the  Gulf  Stream.  This 
is  compared  to  the  predicted  path  of  wave  rays  and  orthogonals  over  the  Gulf  Stream. 
A  number  of  different  possible  current  profiles  for  the  Gulf  Stream  are  examined.  By 
correlating  the  predicted  wave  ray  paths  and  wave  lengths  with  SEASAT  measured 
spectra,  a  best  fit  current  profile  is  proposed. 


Chapter  2 

Theoretical  Background 

2.1      Theory  of  Wave- Current  Interactions 
2.1.1      Governing  Equations 

The  study  of  the  kinematics  of  wave- current  refraction  really  involves  two  problems: 
(1)  the  direction  of  wave  propagation,  and  (2)  the  changes  in  wave  length.  The  wave  or- 
thogonals  (perpendicular  to  wave  crests)  follow  paths  tangent  to  the  wave  number  vector 
at  a  phase  speed  modified  by  current  interactions.  This  is  the  direction  of  propagation 
seen  by  an  external  observer  who  takes  a  'snapshot'  of  the  wave  system.  The  wave  ray 
gives  the  direction  of  energy  propagation  traveling  at  the  group  velocity  and  is  given  by 
the  vector  sum  of  the  intrinsic  group  velocity  and  the  current  velocity  vectors. 

The  assumptions  made  in  this  study  in  applying  the  governing  equations  for  wave- 
current  interactions  are  as  follows: 

1.  The  effects  of  viscosity  and  surface  tension  are  neglected. 

2.  The  fluid  is  assumed  to  be  incompressible  and  its  motion  nearly  irrotational. 

3.  The  currents  are  assumed  to  be  strong  and  vertically  uniform. 

4.  The  spatial  and  temporal  variations  of  the  horizontal  currents  are  assumed  to  be 
slow  relative  to  the  wave  scales;  the  fluid  is  deep  relative  to  the  wave  length  (short- 
wave approximation). 
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CURRENT  NORMAL 
CREST  (NORMAL  TO  ORTHOGONAL) 

Figure  2-1:  Wave  ray  and  orthogonal  relationships 

5.  Reflection  and  diffraction  effects  are  excluded. 

6.  The  wave  amplitude  is  assumed  to  be  small. 

7.  Effects  by  the  surface  curvature  of  the  earth  are  small  enough  to  be  neglected. 

For  a  steady,  inhomogeneous  medium,  the  wave  frequency  is  constant  and  hence  the 
absolute  phase  speed  (c0)  can  be  found  from  the  relative  or  intrinsic  frequency  (u;r), 
wave  number  vector  k,  and  current  vector  U,  i.e., 


a;r  +  k  •  U 
Cr  -f  U  sin  a 


(2.1) 
(2.2) 


where  a  is  the  angle  between  the  wave  orthogonal  (or  the  wave  number  vector  k)  and 

the  current  normal,  as  defined  in  Figure  2-1. 

The  path  of  wave  orthogonals  is  determined  from  integration  of  the  equations: 

dx 
ds 
dy 


ds 
d$ 

ds 


sin# 


e~dt 


[By 


(2.3) 
(2.4) 

(2.5) 
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dx 

du 

dt 

~dk 

dk 

di 

dt 

d 

where  Ca  =  ^  is  the  absolute  phase  speed,  as  seen  by  a  stationary  observer,  9  is  the 
angle  between  the  orthogonal  and  the  x-axis,  and  an  elemental  distance  along  the  arc 
is  given  by  ds  =  cadt.  The  relationship  between  wave  crest  and  current  is  illustrated  in 
Figure  2-1. 

The  wave  ray  path  is  determined  from  integration  of  the  ray  equations  (Kenyon, 
1971): 

(2.6) 
(2.7) 

Here  ^  *s  tne  group  velocity;  gj  describes  the  change  in  wave  number  along  rays. 
Another  form  of  the  ray  equations  (Treloar,  1985)  is: 

(2.8) 
(2.9) 
(2.10) 

where  r  is  now  the  distance  along  the  ray,  and  the  angle  \i  is  defined  by: 

c^  sin  6  +  v 

tan/i  =  -2 (2.11) 

cffrcos0  +  u  v        ' 

where  c^  is  the  relative  group  velocity,  9  is  the  direction  of  the  wave  orthogonal,  and  u 
and  v  are  the  x  and  y  components  of  current  U.  The  relationship  between  the  current 
and  wave  ray  is  also  shown  in  Figure    2-1. 

2.1.2     Analytical  Solutions 

Two  analytical  solutions  are  described  here  to  demonstrate  the  applicability  of  these 
equations.  These  solutions  are  used  to  validate  the  numerical  method  outlined  in  2.1.3. 

Johnson's  Derivation  of  Snell's  Law 

Johnson  (1947)  derived  the  path  of  wave  orthogonals  entering  a  current  from  still  wa- 
ter. This  situation  is  illustrated  in  Figure  2-2.  The  directions  of  wave  crest  propagation 


dx 

dT 

- 

COS/i 

dy 

dT 

= 

sin/i 

dfl 

dr 

= 

,      dtan/x 

CONSTANT  CURRENT 


Figure  2-2:  Wave  orthogonal  entering  a  constant  current  from  still  water 

and  currents  are  related  by  the  following  expression: 

—  =  U  +  -r^r  (2.12) 

sin  a  snap 

The  initial  or  still  water  phase  velocity  is  c,,,  and  the  phase  velocity  relative  to  the 

current  is  c.  The  waves  must  be  continuous  across  the  discontinuity,  so  the  wave  length 

in  still  water  L0  is  related  to  the  wave  length  in  the  current  L: 

L  L0 


sin/} 
Equation  (2.12)  can  be  expressed  as: 

sin/3  =  ■ 


■  Us 


(2.13) 


(2.14) 


Because  in  deep  water  phase  speed  is  related  to  wave  length  by  c2  =  |j,  where  g  is  the 
gravitational  constant,  (2.13)  can  be  rewritten  as: 

(2.15) 


(2.16) 


sm.fi       sina 
A  form  of  Snell's  law  can  be  derived  from  (2.14)  and  (2.15).  Upon  simplifying: 

.  sin  a 

sin/3  = Tr 

P       (1-gsina)2 

The  final  result  (2.16)  describes  the  refraction  of  wave  crests  from  still  water  across  a 

current  discontinuity. 
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Figure  2-3:  Wave  crest  entering  a  current  which  has  a  constant  shear 
Kenyon's  Solution  for  the  Case  of  a  Linear  Shear  Current 

Kenyon  (1971)  used  the  ray  equations  to  derive  an  analytical  solution  to  both  ray 
and  orthogonal  paths  for  the  case  of  a  current  with  a  linear  shear.  The  ray  equations 
hold  under  the  geometrical  optics  approximation  that  wave  amplitude  and  frequency 
vary  slowly  over  distances  of  the  order  of  the  wave  length. 

An  approximate  equation  for  the  radius  of  curvature  of  the  rays  is: 


R-- 


(2.17) 


where  £  =  |(§|  -  ^),  which  is  the  component  of  vorticity  in  the  positive  z  direction. 

The  current  velocity  is  assumed  to  be  in  the  x-direction  and  only  a  function  of  y;  it 

has  a  constant  shear  s: 

sy    ify>0 

0      if  y  <  0 

The  initial  conditions  are  that  the  rays  pass  through  x  =  0,  y  =  0  at  time  t  =  0  with 

initial  angle  <f>0  =  cot-1  -^-.  This  situation  is  shown  in  Figure  2-3. 

The  integration  of  the  ray  equations  (2.6-2.7)  using  the  above  initial  conditions  results 


U{y)-- 


(2.18) 


*y(0      =      KJt  (2.19) 

y(t)     =     2i20(sino0)-JVl-f'2)t -(sino0)-t;  (2.20) 

coto0(sinp0)_*  -  2^7(0,  A")  -  F(eo,K)]}  (2.21) 

where  t'  =  cot  d0-st,  F(Q,K)is  the  elliptic  integral  of  the  first  kind,  cos  0  =  (1-t  2)~t. 
cos0o  =  (sin^0)j,  and  K  =  7T*.  In  (2.20)  and  (2.21)  i^  =  -jf-.  This  is  the  exact 
initial  radius  of  curvature  of  the  rays,  where  va  is  the  initial  wave  frequency. 

Kenyon  also  showed  that  Johnson's  form  of  Snell's  law  can  be  put  in  the  form  x  =  z{y) 
for  the  same  shear  current  by  using  tan  o  =  |f .  This  results  in  a  description  of  the  wave 
orthogonals: 

l(y)  =  'M^We.*)  -  f(6^jt)]  (2-22) 

where  cos  0  =  (sin$0)»(l  +  ^ygf*)"1-  In  Kenyon  (1971),  0  is  misdefined  as  cos2  0  = 

2.1.3     Numerical  Solutions 

Over  the  past  four  decades  numerous  graphical  and  numerical  methods  have  been 
developed  for  refraction  calculations  (see  for  example,  Chao  (1972)  and  Earle  and  Mads  en 
(1987)).  A  somewhat  different  algorithm  has  been  first  described  by  Pararas-Carayannis 
et  aL  (1968)  and  later  by  Abernethy  and  Gilbert  (1975).  Both  algorithms  are  based  on 
triangular  grid  systems  which  assume  that  within  each  triangle  the  velocity  field  is  well 
approximated  by  a  plane.  For  this  condition,  the  ray  path  can  be  solved  analytically 
and  in  particular,  traces  out  the  arc  of  a  circle.  There  are  three  major  advantages  in 
assuming  Vc  constant  within  a  triangular  element  over  the  conventional  approach  with 
a  rectangular  grid.  One,  an  exact  solution  of  (2.5)  is  possible  for  triangles  but  not  for 
rectangular  elements.  Two,  iterative  predictor-corrector  techniques  are  usually  required 
to  advance  the  ray  in  short  steps  across  a  rectangular  grid  element.    Three,  the  ray 
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CENTER  OF  CIRCLE 


Figure  2-4:  Ray  path  along  a  constant  circular  arc  through  a  triangle 


path  calculations  in  triangles  are  much  faster  because  there  is  no  need  to  specify  a  time 
step  nor  a  desired  level  of  accuracy  for  the  iterative  procedure.  Abernethy  and  Gilbert 
(1975)  used  equilateral  triangles  in  their  algorithm.  An  example  of  the  application  of  this 
algorithm  can  be  found  in  Graber  et  al.  (1989),  where  the  variability  of  wave  kinematics 
and  dynamics  due  to  depth  refraction  in  the  North  Sea  is  examined.  Brampton  (1977) 
modified  Abernethy  and  Gilberts'  algorithm  by  using  right  triangles.  These  elements 
allow  variable  grid  spacing  in  the  two  dimensions  and  simplify  the  angular  relationships 
within  the  triangular  grid.  A  typical  triangular  element  can  be  seen  in  Figure  2-4. 

For  this  study  the  algorithm  used  by  Brampton  is  further  modified  to  include  the 
effects  of  currents  in  deep  water.  The  current  field  is  represented  by  discrete  current 
vectors  at  the  vertices  of  triangular  elements  as  seen  in  Figure  2-4.  The  orthogonal 
equations  (2.3),  (2.4),  and  (2.5)  are  applied  to  the  elements.  Using  (2.16)  the  absolute 


wave  velocity  at  a  point  in  a  current  can  be  calculated  from: 

ca  = p- (2.23) 

(l-^cos7)2  V        ' 

where  7  is  the  angle  between  the  wave  orthogonal  and  the  current  shown  in  Figure  2-1. 
By  inserting  a  current  field  with  the  values  specified  at  each  corner  of  the  triangular 
element,  a  new  phase  speed  plane  can  be  calculated  which  includes  the  current  effects. 
The  circular  paths  traced  out  correspond  to  the  wave  orthogonals. 

The  calculation  of  ray  trajectories  is  more  complex,  since  the  orthogonal  path  across 
each  triangle  must  be  computed  first.  The  group  velocity  at  each  corner  is  calculated 
by  taking  the  phase  speed  resulting  from  the  wave-current  interactions.  The  absolute 
group  velocity  is  approximated  by  adding  the  current  components  normal  to  and  along 
the  group  velocity  vector  at  each  point.  In  deep  water  cg  =  \c.  The  entering  angle  of 
the  ray,  /x,  must  also  be  calculated  for  each  element. 

cga     =     [(Cffr  +  tfcos7)2  +  (tfsin7)2]*  (2.24) 

is    =    7  +  tan-*(       Yn1      )  (2-25) 

Kcgr  +  Ucosf/ 

Each  triangle  must  be  analyzed  twice,  first  to  find  7  from  the  orthogonal  path  and  then 
to  trace  the  wave  ray. 

2.2      Examples  of  Solutions 

In  order  to  validate  the  numerical  solution  technique  described  in  the  previous  section, 
a  number  of  scenarios  have  been  posed  in  which  the  numerical  solutions  can  be  compared 
to  analytical  results. 

The  simple  case  first  addressed  by  Johnson  (1947)  involves  the  passage  of  wave  sys- 
tems from  still  water  across  a  discontinuity  into  a  steady,  uniform  current.  In  the  figures 
in  this  section,  the  path  taken  by  wave  orthogonals  are  marked  by  thin  lines  and  the 
wave  rays  are  shown  by  thick  lines.  Table  2.1  compares  results  from  the  analytical  and 
numerical  solutions. 
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10  sec  wave  entering  a  2  m/s  opposing  current. 
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Figure  2-5:  Current  refraction  examples  of  10  second  waves  entering  2  m/s  currents,  (a) 
Following  current,  a  =  30°.  (b)  Opposing  current,  a  =  45°.  Note  that  a  heavy  line 
represents  the  ray  path  and  a  thin  line  the  orthogonal  path. 
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Figure 

Cvirrent 
(2  m/s) 

a 

Analytical 
Orthogonal 

Numerical 
Orthogonal 

Analytical 
Ray 

Numerical 
Ray 

2-5(a) 
2-5(b) 
2-6 

Following 
Opposing 
Following 

30.° 
315.° 
55.° 

34.8° 
323.5° 
Reflected 

35.° 

323.5.° 

Reflected 

44.6° 

338.4° 

Reflected 

44.9° 

337.8° 

Reflected 

Table  2.1:  Comparison  of  analytical  and  numerical  results  for  waves  entering  a  constant 
current. 


A  10  second  wave  which  enters  a  2  m/s  following  current  at  an  angle  greater  than 
53°  from  the  current  normal  will  be  reflected.  This  angle  is  known  as  the  critical  angle. 
This  can  also  be  seen  by  the  solution  generated  by  numerical  method  in  Figure  2-6. 

A  more  complicated  situation  involves  a  wave  system  which  enters  a  current  with  a 
linear  shear.  Figures  2-7  through  2-9  compare  the  outputs  generated  by  Kenyon's  ana- 
lytical solution  with  the  numerical  solution.  It  can  be  seen  that  the  analytical  solutions 
are  almost  identical  to  the  numerical  solutions  to  each  problem,  validating  the  results 
generated  by  the  numerical  model. 

The  numerical  solution  allows  us  to  examine  situations  which  would  involve  very 
complicated  analytical  solutions.  For  instance  Figures  2-10  and  2-ll(a)  show  the  situa- 
tion of  a  ring  shaped  current  loop  with  a  cos2  <f>  profile.  It  can  be  seen  that  waves  may  be 
trapped,  reflected,  or  simply  propagate  through  such  a  current  system  depending  upon 
the  angle  of  incidence. 

Another  interesting  situation  is  the  refraction  of  a  wave  field  propagating  through 
a  large  scale  eddy  such  as  a  warm  core  ring  often  found  near  the  Gulf  Stream.  The 
modeling  of  such  a  warm  core  ring  has  been  described  in  Mapp  et  al.  (1985).  The  flow 
field  used  to  simulate  this  large  eddy  was  a  steady  circular  jet.  The  radial  profile  of 
tangential  velocity  is  a  power  function  joined  to  a  Gaussian.  This  example  can  be  seen 
in  Figure  2- 11(b).  Orthogonals  which  penetrate  into  the  center  of  the  warm  core  ring 
were  refracted,  but  not  greatly.  Orthogonals  which  grazed  the  left  side  of  the  ring  were 
reflected;  they  essentially  entered  a  following  current  at  greater  than  the  critical  angle. 
One  orthogonal  which  entered  the  right  side  of  the  ring  became  trapped  in  an  opposing 
current  and  was  refracted  almost  90°. 


10  sec  wave  entering  a  2  m/s  following  current 
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Figure  2-6:  Example  of  total  reflection  for  a  10  second  wave  entering  a  2  m/s  following 
current,  a  =  55°.  Note  that  a  heavy  line  represents  the  ray  path  and  a  thin  line  the 
orthogonal  path. 
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Figure  2-7:  Refraction  of  11  second  wave  entering  a  following  shear  current  which  varies 
from  0  to  2  m/s,  a  =  60°.  (a)  Solution  generated  using  Kenyon's  analytic  result,  (b) 
Solution  generated  with  numerical  model.  Note  that  a  heavy  line  represents  the  ray  path 
and  a  thin  line  the  orthogonal  path. 


(a) 


1 1  sec  wave  propagating  through  an  opposing  current. 


(b) 


1 1  tec.  wave  propagating  through  an  opposing  current. 


Figure  2-8:  Refraction  of  11  second  wave  entering  an  opposing  shear  current  which  varies 
from  0  to  2  m/s,  a  =  30°.  (a)  Solution  generated  using  Kenyon's  analytic  result,  (b) 
Solution  generated  with  numerical  model.  Note  that  a  heavy  line  represents  the  ray  path 
and  a  thin  line  the  orthogonal  path. 
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propagating  through  an  opposing  current 


(b) 


1 1  sec  wave  propagating  through  an  opposing  current. 


Figure  2-9:  Refraction  of  11  second  wave  entering  an  opposing  shear  current  which  varies 
from  0  to  2  m/s,  a  =  60°.  (a)  Solution  generated  using  Kenyon's  analytic  result,  (b) 
Solution  generated  with  numerical  model.  Note  that  a  heavy  line  represents  the  ray  path 
and  a  thin  line  the  orthogonal  path. 
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20  second  wave  trapped  in  an  opposing  current 


^,>  10  second  wave  reflected  from  a  following  current 
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Figure  2-10:  Examples  of  a  cos2  <f>  current  profile:  (a)  Wave  trapping  (20  second  wave), 
(b)  Wave  reflection  (10  second  wave).  Note  that  a  heavy  line  represents  the  ray  path 
and  a  thin  line  the  orthogonal  path. 


V     \  *  •  • 

-  •                                          t      \  \  •  •    - 

-                                         \  \  '  - 

>  I  * 


(a) 


5  second  wave  propagating  through  an  opposing  current. 
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Figure  2-11:  Examples  of:  (a)  5  second  wave  propagating  through  a  cos2  <j>  current 
profile,  (b)  11.3  second  wave  propagating  through  a  warm  core  ring.  Note  that  a  heavy 
line  represents  the  ray  path  and  a  thin  line  the  orthogonal  path. 


Chapter  3 

SAR  Imaging  of  Ocean  Waves 


3.1      SEASAT 

One  of  the  instruments  deployed  on  the  oceanographic  research  satellite  SEASAT  was 
a  synthetic  aperture  radar.  This  microwave  sensor  artificially  synthesizes  an  aperture  or 
antenna  which  is  hundreds  of  meters  long  in  space.  Therefore,  a  fine  spatial  resolution  of 
about  20  m  can  be  obtained  even  from  altitudes  of  800  km.  The  basic  imaging  response 
of  a  SAR  to  the  ocean  surface  is  explained  by  Bragg  scattering  from  short  gravity  waves. 
In  the  case  of  radars  in  the  frequency  range  of  SEASAT  SAR  (1.27  GHz,  or  the  L- 
Band),  the  Bragg  waves  are  of  the  order  of  20  cm.  The  Bragg  waves,  which  are  formed 
by  local  winds,  are  superimposed  on  underlying  long  (greater  than  fifty  meter)  waves 
which  modulate  the  Bragg  waves.  Thus,  the  Bragg  waves  allow  the  SAR  imaging  of  the 
longer  waves.  This  implies  that  one  requirement  for  SAR  imaging  of  ocean  swell  is  a 
local  wind  field  sufficiently  strong  to  create  these  waves  of  about  20  cm  length. 

One  potential  use  oceanographers  see  for  SAR  is  in  the  collection  of  directional  ocean 
wave  spectra.  Kinsman  (1965)  and  Phillips  (1980)  have  shown  that  the  directional  wave 
energy  spectrum  is  probably  the  most  complete  single  descriptor  of  the  ocean  surface. 
When  the  only  source  of  these  spectra  is  from  local  measurements  taken  from  bouys  or 
ships,  their  use  to  oceanographers  who  seek  to  understand  large  scale  ocean  processes 
is  limited.    Remote  sensing  using  SAR  promises  to  provide  a  valuable  tool  to  better 
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understand  the  ocean  environment. 

The  resolution  obtained  by  S  AR  consists  of  range  resolution  (ground  resolution  per- 
pendicular to  the  satellite  trajectory)  and  azimuth  resolution  (ground  resolution  in  the 
direction  of  the  satellite  trajectory). 

The  range  resolution,  6y,  is  obtained  by  the  pulse  length,  6r,  projected  onto  the 
surface. 

*_-*L_     "L  (3.!) 

2  sine       2  sine  v      ' 

where  c  is  the  velocity  of  light,  r  is  the  pulse  duration,  and  e  is  the  incidence  angle. 

The  azimuthal  resolution  Sx  is  linked  to  the  accuracy  in  determining  the  Doppler 
shift  of  the  backscattered  signal.  A  series  of  pulses  are  transmitted  by  the  radar;  these 
are  sampled  along  the  radar  trajectory,  the  synthesized  aperture.  For  a  stationary  target, 
the  azimuthal  resolution  only  depends  on  the  size  of  the  real  antenna. 

Sm  =  |  (3.2) 

where  D  is  the  real  antenna  length. 

The  imaging  of  waves  by  SAR  is  influenced  by  three  effects.  These  are:  (1)  tilt 
modulation,  (2)  hydrodynamic  interactions,  and  (3)  motion  effects  (velocity  bunching 
and  azimuthal  image  smearing).  These  are  discussed  in  Monaldo  and  Lyzenga  (1986) 
and  in  Graber  and  Kelly  (1988).  Tilt  modulation  is  a  geometric  effect  which  results 
from  the  different  incidence  angles  of  Bragg  scattering  waves  seen  by  the  radar  as  they 
propagate  over  longer  waves.  The  hydrodynamic  effect  is  characterized  by  modulations 
in  power  caused  by  changes  in  the  spectral  density  at  the  Bragg  wave  number  across  the 
profile  of  the  long  waves.  Hydrodynamic  modulations  result  from  long  wave-short  wave 
interactions  which  increase  the  ripple  amplitude  near  swell  crests  and  decrease  the  ripple 
height  near  swell  troughs.  They  also  come  from  the  distortion  of  the  airflow  above  the 
sea  surface  by  swell  waves.  These  effects  are  environmentally  dependent  and  generally 
not  well  understood.  Velocity  bunching  causes  image  degradation  due  to  the  movement 
of  the  sea  surface  resulting  from  the  orbital  velocities  of  long  waves.  The  rising  face  of 
the  wave  is  shifted  on  the  image  in  the  direction  of  flight,  while  the  falling  face  is  shifted 
opposite  the  flight  direction.   For  waves  traveling  in  range,  this  has  no  imaging  effect. 
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However,  waves  propagating  in  the  azimuth  direction  are  affected  significantly.  Similar 
effects  caused  by  waves  shorter  than  the  azimuthal  resolution  dimension  cause  Doppler 
smearing. 

The  combination  of  these  two  effects  creates  a  limitation  faced  by  all  SARs  viewing 
the  ocean:  a  degradation  of  the  wave  energy  spectrum  in  the  azimuth  direction.  The 
degradation  is  particularly  pronounced  in  high  flying  SARs  such  as  the  one  deployed 
on  SEASAT.  Azimuthal  smearing  limits  the  usefulness  of  SEASAT  data  in  studying 
the  refraction  of  wave  systems  across  currents.  As  is  discussed  in  section  4.3,  the  wave 
system  examined  in  this  case  study  is  refracted  from  range  into  azimuth.  This  fact  made 
it  difficult  to  track  the  wave  system  across  the  Gulf  Stream. 

3.2      Orbit  1339  Images 

The  synthetic  aperture  radar  deployed  on  SEASAT  had  a  lifetime  of  only  about  three 
months.  There  are  a  limited  number  of  satellite  passes  which  had  interesting  wind  and 
wave  conditions  along  with  independent  measures  of  these  parameters.  Pass  1339  of 
September  28,  1978  was  one  such  pass. 

Beal  et  al.  (1983)  and  (1986)  reported  on  the  processing  of  the  SAR  images  and  the 
oceanographic  conditions  along  the  entire  pass.  The  present  study  focuses  on  the  portion 
of  the  pass  in  the  vicinity  of  the  Gulf  Stream.  Figures  3-1  and  3-2  which  appeared  in 
Beal  et  al.  (1986)  show  the  geographic  orientation  of  pass  1339. 

Digitally  processed  SAR  data  from  MacDonald  Dettwiler  and  Associates  consisting 
of  40  i  40  km  images  were  used  for  the  analysis  here.  A  description  of  the  processing  of 
these  images  is  included  in  Chapter  4. 


Figure  3-1:  SEAS  AT  pass  1339  over  the  North  Atlantic  (from  Beal  et  al.  : 
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Figure  3-2:  SEASAT  pass  1339  over  the  Gulf  Stream  (from  Beal  et  al.  1986). 
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Chapter  4 

Case  Study:   Gulf  Stream 

4.1      Environmental  and  Meteorological  Conditions  Dur- 
ing Pass  1339 

SEASAT  pass  1339  crossed  the  Gulf  Stream  on  September  28,  1978  at  1521  GMT. 
Two  dominant  wave  systems  existed  throughout  most  of  the  900  km  pass.  The  primary 
system  consisted  of  200  m  swell  which  propagated  from  southeast  to  northwest;  a  sec- 
ondary system  of  100  m  waves  was  more  variable  in  direction,  but  generally  traveled 
between  a  westerly  to  southwesterly  direction. 

Beal  et  al.  (1986)  demonstrated  that  the  primary  system  originated  in  Tropical 
Depression  12  on  September  24.  The  maximum  winds  in  the  vicinity  of  the  depression 
were  estimated  at  20  m/s  and  generated  waves  which  traveled  towards  the  west.  This 
wave  source  can  be  seen  in  Figure  4-1. 

A  smaller  storm  closer  to  the  satellite  pass  developed  on  26  and  27  September.  The 
peak  winds  were  estimated  at  15  m/s.  This  storm  can  also  be  seen  on  Figure  4-1  and  is 
the  source  of  the  secondary  100  m  wave  system. 

In  the  immediate  vicinity  of  the  satellite  pass,  a  weak  trough  dominates  the  local 
wind  field  causing  a  counter  clockwise  wind  movement  over  the  900  km  length  of  the 
satellite  pass.  The  wind  is  southerly  with  a  speed  of  about  5  m/s  over  the  Gulf  Stream. 
The  wind  vectors  are  shown  in  Figure  3-1. 


Figure  4-1:  Sources  of  wave  systems  observed  by  SEASAT  SAR  on  pass  1339  (from  Beal 
et  al.  1986). 
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4.2      SAR  Data  Analysis 

The  available  set  of  SAR  data  from  pass  1339  consists  of  nineteen  40  x  40  km  images 
digitally  processed  by  MacDonald  Dettwiler  and  Associates  (Vanvouver,  Canada).  Beal 
et  al.  (1983)  and  (1986)  looked  at  the  entire  pass  to  follow  the  propagation  of  the  200 
m  swell  system.  In  the  current  study,  only  the  two  40  km  images  which  span  the  Gulf 
Stream  are  re-examined  in  a  detailed  analysis  of  the  wave-current  interactions,  as  the 
200  m  swell  system  enters  the  flow  of  the  Gulf  Stream  at  an  oblique  angle  and  undergoes 
refraction. 

From  each  of  the  two  images,  eleven  frames  of  size  512  x  512  pixels  were  defined. 
These  provided  overlapping  images  of  6.4  x  6.4  km  on  a  side.  Each  of  the  twenty-two 
frames  thus  defined  overlaps  the  next  by  about  3  km  with  the  exception  of  the  first  and 
last  in  each  of  the  large  images. 

The  data  from  each  frame  is  normalized  by  subtracting  the  mean  value  of  the  frame 
from  each  pixel  value  within  the  frame  and  then  dividing  by  the  mean.  A  two  dimensional 
fast  Fourier  transform  is  performed  on  each  of  the  resulting  normalized  data  frames.  The 
complex  transform  is  multiplied  by  its  complex  conjugate  and  oriented  with  the  zero  wave 
number  value  in  the  center.  The  resulting  transform  contains  information  out  to  §f  ^- 
The  image  is  then  trimmed  to  272  x  272  pixels.  An  example  of  a  raw  Fourier  transform 
is  shown  in  Figure  4-2. 

The  transforms  are  spectrally  smoothed  using  a  discrete  15  x  15  pixel  Gaussian 
function.  Following  Beal  et  al.  (1986),  the  smoothed  spectrum  St  is  related  to  the  raw 
spectrum  5U  by: 

«.(*..*,)- +7     +7  -S2 <") 

Lm=-7Lj=-7€        " 

The  effect  of  the  smoothing  kernel  (4.1)  on  the  raw  spectrum  is  also  shown  in  Figure  4- 
2.  In  this  and  following  figures  showing  contour  plots,  adaptive  contour  intervals  were 
used  to  highlight  energy  peaks  of  interest.  The  interior  of  the  256  x  256  pixel  smoothed 
region  represent  wave  lengths  greater  than  50  m. 

The  spectra  could  be  further  smoothed  by  along-track  averaging  of  several  adjacent 
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(a) 


Unsmoothed  Spectrum 


(b) 


Smoothed  Spectrum 
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Figure  4-2:  (a)  Raw  FFT  of  image  20  which  is  located  just  after  the  Gulf  Stream,  (b) 
Smoothed  FFT  (same  image  as  (a)).  Satellite  travel  (azimuth  direction)  is  along  the 
x-axis.  ni 


spectra,  but  this  is  inappropriate  over  a  region  of  rapidly  changing  conditions  such  as 
the  Gulf  Stream. 

4.3      Results 

Beal  et  al.  (1986)  processed  six  6.4  x  6.4  km  spectra  from  each  40  x  40  km  SAR 
image.  As  they  followed  what  appeared  to  be  the  dominant  200  m  wave  system  across 
the  Gulf  Stream,  they  observed  a  wave  shortening  resulting  in  a  wave  number  increase 
four  times  that  which  their  model  predicted.  The  wave  number  vector  seemed  to  turn 
into  the  current,  contrary  to  the  behavior  predicted  by  Snell's  law.  Beal  et  al.  concluded 
that  it  is  possible  that  the  system  they  traced  may  actually  be  a  secondary  system  inside 
the  Gulf  Stream.  Figure  4-3  shows  the  wave  length  and  angle  of  the  system  which  Beal  et 
al.  traced  across  the  Gulf  Stream.  The  wave  length  behaved  approximately  as  expected 
for  a  wave  which  was  propagating  into  an  opposing  current  at  an  oblique  angle,  but  the 
direction  turned  in  the  wrong  sense. 

In  this  study  a  closer  spacing  of  wave  spectra  is  used  to  allow  trends  to  be  better 
seen  in  the  rapidly  changing  environment  of  the  Gulf  Stream.  It  can  be  seen  that  merely 
chosing  the  highest  energy  spectral  peak  does  not  ensure  that  the  correct  system  is  always 
observed  during  the  wave  system's  propagation  through  the  Gulf  Stream.  Figure  4-4 
show  the  geographical  orientation  of  the  SEASAT  SAR  images  from  pass  1339  which 
were  processed  into  the  twenty-two  spectra.  Figures  4-5  through  4-10  display  these 
spectra.  The  energy  peak  of  the  200  m  wave  system  is  marked  on  each  figure.  It  can  be 
seen  that  azimuthal  smearing  makes  the  determination  of  this  peak  difficult  inside  the 
Gulf  Stream-  The  wave  lengths  and  direction  of  these  peaks  are  plotted  on  Figure  4- 
11.  The  shortening  and  refraction  of  the  wave  system  toward  the  North  inside  the  Gulf 
Stream  is  visible.  These  results  suggest  that  Beal  et  al.  did  apparently  trace  the  wrong 
system  through  the  Gulf  Stream. 

It  can  be  seen  from  Figure  4-11  that  the  wave  length  of  the  wave  system  across  the 
Gulf  Stream  is  shortened  from  200  m  to  about  170  m.  At  a  40°  incidence  angle  between 
the  wave  front  and  the  Gulf  Stream  normal,  a  200  m  wave  should  be  shortened  to  175  m 
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Figure  4-3:  (a)  Wave  length  and  (b)  wave  angle  variation  across  the  Gulf  Stream  (from 
Beal  et  al.  1986). 
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Figure  4-4:  Geographical  orientation  of  SEAS  AT  pass  1339  over  the  Gulf  Stream.  Num- 
bers at  right  of  track  indicate  high- resolution  SAR  image  sequence  1  through  22,  and 
numbers  at  left  refer  to  the  two  larger  40  x  40  km  images. 
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Figure  4-5:  SAR  wave  spectra:  images  1  through  4. 
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Figure  4-6:  SAB.  wave  spectra:  images  5  through  8.  The  Gulf  Stream  edge  is  located 
between  images  7  and  8.  Note  that  the  energy  peak  corresponding  to  the  wave  system 
of  interest  is  not  evident  in  image  5. 
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Figure  4-7:  SAR  wave  spectra:  images  9  through  12.  Note  a  6.7  km  gap  occurs  between 
spectra  11  and  12. 
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Figure  4-8:    SAB.  wave  spectra:    images  13  through  16.    The  North  Wall  of  the  Gulf 
Stream  is  located  between  images  15  and  16. 
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Figure  4-9:  SAR  wave  spectra:  images  17  through  20. 
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Figure  4-10:  SAR.  wave  spectra:  images  21  and  22. 
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Angle  of  wave  system  from  SAR 


Distance  along  x-axis  (m) 


Figure  4-11:  Wave  length  and  direction  of  wave  front  crossing  the  Gulf  Stream.  The 
error  bars  represent  the  range  of  possible  values  along  the  ridges  seen  in  the  spectra  in 
the  Gulf  Stream.  The  x-axis  is  oriented  along  the  east-west  direction;  the  waves  are 
propagating  toward  the  coast. 
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in  the  presence  of  a  2.0  m/s  current  and  to  165  m/s  in  the  presence  of  a  2.5  m/s  current. 
This  is  in  the  range  of  current  velocities  often  observed  in  the  Gulf  Stream  (Williams 
et  al.  1974).  Within  the  Gulf  Stream,  azimuthal  smearing  results  in  the  energy  center 
appearing  as  a  ridge  rather  than  as  a  clear  peak.  The  range  of  possible  values  along 
this  ridge  are  shown  by  the  error  bars  within  the  Gulf  Stream  on  Figure  4-11.  The  wave 
lengths  are  more  variable  than  the  angles,  because  the  smearing  in  azimuth  degrades  the 
wave  lengths  more  than  the  angles. 

Figures  4-12  and  4-13  show  the  SAR  images  reduced  in  resolution.  These  allow  us 
to  visually  examine  the  southern  and  northern  boundaries  of  the  Gulf  Stream.  The 
southern  boundary  seen  in  Figure  4-12  falls  at  an  angle  of  about  7°  or  8°  relative  to 
the  satellite  trajectory.  Since  the  satellite  pass  was  23°  from  North,  this  indicates  that 
the  current  was  flowing  at  about  75°  from  North  at  the  southern  boundary  of  the  Gulf 
Stream.  The  northern  boundary  is  very  distinct  and  can  be  measured  to  rise  at  an  angle 
of  7°.  This  would  correspond  to  a  current  direction  of  60°  at  the  North  Wall  of  the  Gulf 
Stream.  These  measurements  indicate  that  the  stream  is  widening  at  the  point  at  which 
the  satellite  passes  over. 

The  data  shown  in  Figure  4-11  imply  that  the  refraction  effect  is  strongest  near  the 
northern  boundary  of  the  Gulf  Stream.  This  fact  would  suggest  an  asymmetric  current 
with  a  velocity  peak  shifted  towards  the  North  Wall.  This  might  also  be  a  reason  why 
the  boundary  is  more  sharply  denned  in  Figure  4-13  than  in  Figure  4-12. 

In  order  to  use  the  SAR  data  to  estimate  the  current  field  at  the  point  at  which  the 
satellite  passed  over  the  Gulf  Stream,  a  number  of  possible  currents  were  examined  with 
the  refraction  model.  The  width  of  the  Gulf  Stream  current  is  estimated  to  be  33.6  km. 
This  is  based  on  measurements  from  the  40  x  40  km  images  as  well  as  from  the  eight 
spectra  between  the  start  and  end  of  the  Gulf  Stream  overpass.  The  spectra  have  an 
average  spacing  of  3.36  km,  and  there  is  another  gap  of  that  size  between  the  40  x  40  km 
images.  The  satellite  pass  direction  is  nearly  normal  to  the  Gulf  Stream,  so  a  total  of 
ten  3.36  km  spacings  results  in  a  current  approximately  34  km  wide.  The  grid  spacing 
used  in  the  numerical  model  is  51  x  51  grid  points  each  1.68  km  apart  (i.e.,  the  model 
grid  spacing  is  half  the  size  of  the  spectra  spacing).  The  assumption  is  also  made  that 
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Figure  4-12:  First  40  i  40  km  SAR  image  showing  the  location  of  the  southern  edge  of 
the  Gulf  Stream  (arrow). 
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Figure  4-13:  Second  40  x  40  km  SAR  image.  The  northern  boundary  (North  Wall)  of 
the  Gulf  Stream  is  clearly  visible  in  the  image  (arrow). 
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the  current  direction  is  constant,  so  eddies  and  lateral  variations  within  the  Gulf  Stream 
are  discounted. 

Current  directions  were  proposed  to  vary  from  60°  to  75°.  The  current  profiles  ex- 
amined vary  with  a  cosine- squared  dependence  (symmetic  jet),  a  skewed  cosine-squared 
dependence  with  the  peak  closer  to  the  northern  boundary  of  the  Gulf  Stream  (asym- 
metric jet),  and  a  top  hat  shape  with  tapered  edges  (top  hat).  Cross-sections  of  these 
profiles  can  be  seen  in  Figure  4-14.  Peak  current  speeds  (Umax)  within  these  profiles  of 
1.5,  2.0,  2.25  and  2.5  m/s  were  examined.  Each  of  these  current  directions,  profiles,  and 
speeds  were  input  into  the  numerical  model  for  a  200  m  wave  with  an  initial  direction 
of  295°.  Note  that  the  average  angle  of  the  center  of  the  spectra  before  the  Gulf  Stream 
was  295°,  and  the  average  wave  length  was  200  m  (11.3  second  period).  Figures  4-15 
through  4-17  show  the  orthogonals  and  rays  being  refracted  across  six  model  currents. 

The  wave  rays  show  the  propagation  of  wave  energy,  so  they  are  compared  to  the 
spectra  processed  from  the  SAR  data.  Figures  4-18  and  4-19  show  ray  angles  and  wave 
lengths  from  the  SAR  data  (+)  superimposed  on  the  model  output.  Both  the  model  and 
the  data  are  oriented  along  the  x-axis.  For  each  data  point,  the  corresponding  model 
values  for  angle  and  wave  length  are  found  from  simple  interpolation.  From  these  two 
series  of  wave  lengths  and  angles,  two  sets  of  vectors  are  defined  in  the  direction  of  the 
wave  ray  and  with  the  magnitude  of  the  wave  length.  The  SAR  data  are  compared  to  the 
model  output  by  finding  the  phase  angle  of  the  complex  correlation  coefficient  between 
the  two  vector  series.  This  technique  is  described  in  Kundu  (1975).  The  vectors  are 
defined  as: 

u(t)  =  u{t)  +  iv{t)  (4.2) 

and  are  divided  into  east,  u(t),  and  north,  v(t),  components.  The  complex  correlation 
coefficient  between  the  model  vector  series  u>m(t)  and  the  SAR  data  vector  series  a;j(t) 

<aM')> (4.3) 

<  u)m(t)u;d(t)  >2<  u>Z{t)um{t)  >* 

where  v'{t)  is  the  complex  conjugate.  The  correlation  can  also  be  expressed  in  terms  of 


east-north  components,  i.e., 


<  "mtij  +  vmvd  >  .  <urnvd-udvm> 

—r  +  I j r  (4.4) 


<  u^  +  t/^  >  a  <  ud  +  t£  >5         <  u^  +  1&  >  2  <  ud  +  vd  >  * 
and  the  phase  angle  is  : 

-1   <  UmVd  ~  UdVm  > 


(4.5) 


Note  that  <  ...  >  indicates  the  average  value  of  the  included  series.  These  complex 
correlations  are  listed  in  Tables  4.1,  4.2,  and  4.3.  All  magnitudes  of  correlation  are  high 
because  the  vectors  are  all  pointing  close  to  the  same  direction. 

To  provide  another  measure  of  how  well  the  numerical  model  output  fits  the  SAR 
data,  the  wave  length  and  angle  are  examined  individually  by  calculating  residuals.  The 
residual  variance  in  wave  length  and  angle  are  calculated  from  the  deviation  of  the  model 
results  from  the  SAR  data: 

wave  length  residual  =  —  ^(i„,|  -  m^j)2  (4-6) 

angle  residual  =  —  ^{xang  -  mang)2  (4.7) 

(4.8) 

where  N  is  the  number  of  data  points  xwi  and  xang  are  the  SAR  measured  wave  lengths 
and  angles,  and  m^i  and  mang  are  the  numerical  model  generated  wave  lengths  and 
angles. 

The  complex  correlations  and  residuals  are  examined  together  to  find  the  best  overall 
fit  between  the  SAR  data  and  the  model  output  in  order  to  determine  the  current 
characteristics  which  could  explain  the  observed  wave  refraction. 

As  can  be  seen  from  Tables  4.1,  4.2,  and  4.3,  the  best  overall  correlation  between  the 
numerical  refraction  model  and  the  SAR  data  is  found  when  the  current  is  represented 
by  a  top  hat  profile,  a  velocity  of  2.0  m/s,  and  a  direction  75°.  In  general,  the  correlation 
improves  as  the  profile  is  changed  from  a  symmetric  jet,  to  an  asymmetric  jet,  to  a  top 
hat  shape.  This  indicates  that  most  of  the  refraction  occurs  as  the  wave  enters  and  leaves 
the  Gulf  Stream,  so  the  current  profile  within  the  Gulf  Stream  at  the  location  of  the 
satellite  pass  is  close  to  uniform,  and  the  largest  velocity  gradients  are  near  the  current 
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Figure  4-14:  Model  current  profiles:  (a)  Symmetric  jet.   (b)  Asymmetric  jet  .   (c)  Top 
hat. 
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(a)  10 


x  104  1.5  m/s  symmetric  jet  current  at  65  deg. 


(b) 


2.5  m/s  symmetric  jet  current  at  75  deg. 


x-axis  (m) 


Figure  4-15:  Refraction  of  wave  rays  and  orthogonals  across  symmetric  jet  current  pro- 
files, (a)  Umax  is  1.5  m/s  at  65°.  (b)  Umax  is  2.5  m/s  at  75°.  Note  that  a  heavy  line 
represents  the  ray  path  and  a  thin  line  the  orthogonal  path. 
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(a) 


xlO*  1.5  m/s  asymmetric  jet  current  at  60  deg 


(b) 


xlO4  2.5  m/s  asymmetric  jet  current  at  75  deg 


x-axis  (m) 


Figure  4-16:  Refraction  of  wave  rays  and  orthogonals  across  asymmetric  jet  current 
profiles,  (a)  Umax  is  1.5  m/s  at  60°.  (b)  Umax  is  2.5  m/s  at  75°.  Note  that  a  heavy  line 
represents  the  ray  path  and  a  thin  line  the  orthogonal  path. 
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(a) 


1.5  m/s  top  hat 


at  75  deg 


(b) 


2.0  m/s  top  hat  current  at  60  deg 


Figure  4-17:  Refraction  of  wave  rays  and  orthogonals  across  top  hat  current  profiles,  (a) 
Umax  is  1.5  m/s  at  75°.  (b)  Umax  is  2.0  m/s  at  60°.  Note  that  a  heavy  line  represents 
the  ray  path  and  a  thin  line  the  orthogonal  path. 
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2.0  m/s  asymmetric  jet  coircnt  at  75  deg 


Distance  along  x  axis  (m) 


Figure  4-18:  (a)  Symmetric  jet  profile,  (b)  Asymmetric  jet  profile.  SAR  data  (+) 
superimposed  on  model  output;  the  waves  are  propagating  toward  the  coast.  The  x-axis 
corresponds  to  the  east- west  direction. 
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(a)  330 


2.0  m/s  top  hat  current  at  75  deg 


(b) 


Z25  m/s  top  hat  current  at  70  deg 


Distance  along  x  axis  (m) 


Figure  4-19:  (a)  Top  hat  profile  with  2  m/s  current,  (b)  Top  hat  profile  with  2.25 
m/s  current.  SAR  data  (+)  superimposed  on  model  output;  the  waves  are  propagating 
toward  the  coast.  The  x-axis  corresponds  to  the  east- west  direction. 
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Umax 

Current 

Magnitude  of 

Phase 

Angle 

Wave  Length 

Direction 

Correlation 

Angle 

Residual  deg2 

Residual  m2 

1.5 

60° 

.9908 

4.0180° 

89.8 

136.3 

1.5 

65° 

.9916 

3.6014° 

82.1 

125.1 

1.5 

70° 

.9923 

2.9953° 

74.1 

118.2 

1.5 

75° 

.9932 

2.7280° 

67.4 

110.9 

2.0 

60° 

.9922 

3.0911° 

74.6 

127.1 

2.0 

65° 

.9926 

2.3074° 

67.8 

121.8 

2.0 

70° 

.9937 

1.8731° 

60.7 

111.7 

2.0 

75° 

.9952 

1.7376° 

50.5 

113.6 

2.25 

60° 

.9928 

3.0420° 

70.1 

128.8 

2.25 

65° 

.9932 

2.1865° 

63.1 

127.0 

2.25 

70° 

.9943 

1.7043° 

56.7 

114.6 

2.25 

75° 

.9950 

1.1597° 

50.7 

123.8 

2.5 

60° 

.9919 

2.3578° 

73.0 

132.6 

2.5 

65° 

.9928 

1.4241° 

63.0 

133.4 

2.5 

70° 

.9942 

1.3091° 

56.4 

124.8 

2.5 

75° 

.9956 

1.4134° 

48.7 

127.4 

Table  4.1:    Correlations  of  SAR  data  vs  numerical  refraction  model  for  symmetric  jet 
current  profile. 


UmaX 

Current 

Magnitude  of 

Phase 

Angle 

Wave  Length 

Direction 

Correlation 

Angle 

Residual  deg2 

Residual  m2 

1.5 

60° 

.9908 

4.0341° 

89.6 

134.3 

1.5 

65° 

.9918 

3.4259° 

79.7 

121.3 

1.5 

70° 

.9928 

3.0297° 

71.5 

108.9 

1.5 

75° 

.9933 

2.3724° 

64.3 

110.8 

2.0 

60° 

.9927 

3.0574° 

70.4 

123.2 

2.0 

65° 

.9928 

2.0540° 

65.7 

119.2 

2.0 

70° 

.9939 

1.5348° 

57.8 

112.9 

2.0 

75° 

.9958 

1.5629° 

46.8 

106.1 

2.25 

60° 

.9936 

3.1710° 

66.0 

122.1 

2.25 

65° 

.9942 

2.3687° 

58.5 

117.2 

2.25 

70° 

.9953 

1.7400° 

50.2 

110.8 

2.25 

75° 

.9955 

0.6293° 

46.7 

118.3 

2.5 

60° 

.9929 

2.7747° 

68.2 

126.0 

2.5 

65° 

.9938 

1.4798° 

58.0 

120.0 

2.5 

70° 

.9951 

1.6270° 

51.2 

122.5 

2.5 

75° 

.9962 

2.4065° 

48.4 

122.0 

Table  4.2:   Correlations  of  SAR  data  vs  numerical  refraction  model  for  asymmetric  jet 
current  profile. 


Umax 

Current 

Magnitude  of 

Phase 

Angle 

Wave  Length 

Direction 

Correlation 

Angle 

Residual  deg2 

Residual  m2 

1.5 

60° 

.9921 

3.1317° 

75.7 

112.0 

1.5 

65° 

.9940 

2.4323° 

61.7 

95.5 

1.5 

70° 

.9941 

2.1287° 

58.2 

92.8 

1.5 

75° 

.9953 

1.0883° 

49.3 

96.0 

2.0 

60° 

.9930 

1.7350° 

64.3 

108.3 

2.0 

65° 

.9945 

0.9586° 

54.8 

99.7 

2.0 

70° 

.9936 

0.5190° 

48.4 

94.1 

2.0 

75° 

.9963 

0.1357° 

45.1 

98.9 

2.25 

60° 

.9937 

2.0353° 

61.7 

110.2 

2.25 

65° 

.9926 

0.0347° 

66.2 

114.0 

2.25 

70° 

.9965 

1.3584° 

46.4 

106.8 

2.25 

75° 

.9972 

0.5551° 

41.5 

111.2 

2.5 

60° 

.9929 

2.5981° 

70.0 

120.3 

2.5 

65° 

.9943 

2.5297° 

63.4 

119.1 

2.5 

70° 

.9948 

1.8436° 

60.6 

136.0 

2.5 

75° 

.9955 

-.1881° 

54.5 

126.2 

Table  4.3:  Correlations  of  SAR  data  vs  numerical  refraction  model  for  top  hat  current 
profile. 


edges.  The  correlation  improvement  from  the  symmetric  jet  to  the  asymmetic  jet  shape 
indicates  that  the  velocity  may  be  slightly  skewed  towards  the  North  Wall. 

The  current  direction  which  provides  the  best  fit  with  the  SAR  data  seems  to  be 
between  70°  and  75°.  The  southern  boundary  of  the  Gulf  Stream  is  oriented  at  about 
75°,  so  this  value  is  reasonable.  The  Umax  which  provides  the  best  correlation  is  2.25 
m/s  for  the  symmetric  jet  current  profile,  2.5  m/s  for  the  asymmetric  jet  current  profile, 
and  2.0  m/s  for  the  top  hat  profile. 

There  is  a  fair  amount  of  data  scatter  both  before  and  within  the  Gulf  Stream.  This 
makes  it  difficult  to  assign  a  specific  meaning  to  a  pair  of  data  points  which  digress  from 
the  trend  followed  by  the  rest.  However,  both  the  wave  length  and  angle  of  propagation 
show  two  data  points  within  the  Gulf  Stream  which  differ  significantly  from  a  smooth 
curve.  This  could  be  indicative  of  an  eddy  formation  or  lateral  variability  within  the 
Gulf  Stream  at  this  point;  the  length  scale  is  between  6  km  and  12  km. 


Chapter  5 

Summary  and  Conclusions 


Ocean  waves  which,  interact  with  currents  in  deep  water  undergo  refraction,  and 
characteristics  of  the  waves  like  length  and  height  are  affected  by  the  interaction.  This 
type  of  refraction  can  be  observed,  as  swell  propagates  through  major  current  systems 
like  the  Gulf  Stream  or  the  Antrarctic  Circumpolar  Current.  A  model  which  examines 
wave-current  kinematics  can  be  used  in  conjunction  with  remote  sensing  measurements 
of  wave  propagation  to  infer  the  current  system  which  caused  the  measured  refraction. 

The  study  of  the  kinematics  of  wave-current  refraction  really  involves  two  problems: 
(1)  the  direction  of  wave  propagation,  and  (2)  the  changes  in  wave  length.  The  wave  or- 
thogonals  (perpendicular  to  wave  crests)  follow  paths  tangent  to  the  wave  number  vector 
at  a  phase  speed  modified  by  current  interactions.  This  is  the  direction  of  propagation 
seen  by  an  external  observer  who  takes  a  'snapshot'  of  the  wave  system.  The  wave  ray 
gives  the  direction  of  energy  propagation  traveling  at  the  group  velocity  and  is  given  by 
the  vector  sum  of  the  intrinsic  group  velocity  and  the  current  velocity. 

A  numerical  model  is  developed  which  predicts  the  direction  of  propagation  and  wave 
length,  as  a  wave  ray  or  orthogonal  is  refracted  by  a  current.  This  model  is  validated  by 
comparing  its  results  to  analytical  solutions  for  cases  of  simple  refraction  problems. 

SEASAT  SAR  data  from  a  pass  over  the  Gulf  Stream  are  analyzed  to  allow  wave 
energy  propagation  to  be  traced  across  the  Gulf  Stream.  The  digitally  processed  SAR 
images  are  sectioned  into  overlapping  512  x  512  pixel  frames.  These  frames  are  normal- 


ized,  and  a  Fast  Fourier  Transform  is  performed  on  each  one.  The  magnitude  of  the 
transform  is  found,  and  it  is  oriented  with  the  zero  wave  number  at  the  center.  The 
resulting  raw  spectra  are  smoothed  and  trimmed  to  obtain  energy  spectra  which  show 
waves  of  length  greater  than  50  m.  These  spectra  contain  energy  peaks  which  corre- 
spond to  a  200  m  wave  system  which  traveled  from  southeast  to  northwest  across  the 
Gulf  Stream.  The  peaks  are  distinctly  visible  before  and  after  the  Gulf  Stream,  but 
within  the  Gulf  Stream  the  waves  are  refracted  into  azimuth,  and  the  effect  of  azimuth 
or  Doppler  smearing  makes  it  much  more  difficult  to  trace  the  wave  energy  propagation. 

The  general  position  and  width  of  the  Gulf  Stream  below  the  satellite  pass  is  clearly 
marked  in  the  SAR  images.  Three  possible  current  profiles  are  proposed  to  model  the 
Gulf  Stream:  a  symmetric  jet,  an  assymetric  jet,  and  a  top  hat.  The  current  direction  is 
modeled  by  directions  of  flow  which  vary  from  60°  to  75°,  and  the  current  velocity,  Umax, 
is  varied  from  1.5  to  2.5  m/s.  These  current  parameters  are  input  into  the  numerical 
model  as  a  current  field.  Using  the  model,  the  kinematics  of  a  200  m  wave  passing 
through  the  Gulf  Stream  are  examined.  The  model  output  is  compared  to  the  SAR  data 
by  colocating  the  model  with  the  data  by  representing  the  wave  length  and  direction  of 
propagation  as  sets  of  vectors.  From  this  the  magnitude  and  phase  angle  of  the  complex 
correlation  between  the  SAR  data  and  the  model  are  calculated.  The  wave  length  and 
angle  residuals  are  also  calculated  to  provide  an  additional  measure  for  comparison. 

The  best  correlations  are  found  for  Gulf  Stream  velocities  between  2.0  and  2.5  m/s 
at  70°  to  75°.  The  best  overall  correlation  between  the  model  and  SAR  data  is  obtained 
when  the  Gulf  Stream  is  modeled  by  a  top  hat  profile,  a  velocity  of  2.0  m/s,  and  a 
direction  of  75°.  75°  agrees  with  the  direction  visually  observed  from  the  SAR  images, 
and  the  direction  and  speeds  are  close  to  the  Coast  Guard  estimates  for  the  Gulf  Stream 
on  this  day.  The  current  profiles  used  in  the  model  do  not  take  into  account  a  possible 
widening  of  the  Gulf  Stream  at  the  position  of  the  satellite  overpass.  There  is  a  great 
deal  of  scatter  in  the  SAR  data  both  before  and  in  the  Gulf  Stream,  so  it  is  difficult  to 
correlate  every  point  with  specific  current  behavior,  but  the  increase  in  wave  length  and 
change  in  angle  in  the  center  of  the  Gulf  Stream  seem  to  indicate  that  there  may  be  a 
non-uniform  feature  such  as  an  eddy  or  lateral  variation  forming  within  the  current. 
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SAR  data  can  be  used  to  trace  wave  systems  which  pass  through  a  current  like  the 
Gulf  Stream.  The  tomographic  employment  of  SAR  data  in  conjunction  with  a  numerical 
model  for  wave-current  refraction  can  be  used  to  infer  the  current  which  caused  the  wave 
refraction.  In  order  for  the  approach  described  above  to  be  practically  useful,  a  uniform 
and  stable  long  wave  system  must  exist.  This  tomographic  application  requires  a  wave 
system  which  stays  in  the  range  direction  with  regard  to  the  satellite,  or  which  is  long 
enough  to  still  be  measured  in  azimuth.  The  large  scale  current  behavior  is  quite  evident 
in  the  SAR  data,  but  small  scale  eddies  or  lateral  variations  are  harder  to  measure, 
because  they  may  have  smaller  length  scales  than  the  scales  of  the  frames  processed 
here. 

Lower  flying  SARs  which  overcome  some  of  the  inherent  problems  of  azimuth  smear- 
ing are  superior  to  the  SAR  which  flew  in  SEASAT  for  this  type  of  tomographic  ap- 
plication. Two  Shuttle  missions  deployed  SARs,  the  SIR-A  and  SLR-B  flights;  they 
provided  relatively  short  sets  of  SAR  data  from  lower  altitudes  than  SEASAT.  Even 
SARs  deployed  in  high  flying  aircraft  can  provide  wide  geographic  coverage  of  current 
systems. 

SAR  and  numerical  models  of  wave-current  refraction  offer  promising  new  tools  to 
the  oceanographer,  as  he  seeks  to  better  measure  and  understand  the  ocean  environment. 
When  they  are  used  together,  the  structure  of  current  systems  can  be  measured  from 
space.  This  type  of  tomography  is  limited  by  azimuth  smearing,  as  a  wave  system  is 
refracted  out  of  the  range  and  into  the  azimuth  direction.  Lower  flying  SARs  which  will 
overcome  azimuth  smearing  and  also  provide  better  small  scale  resolution  are  needed  to 
continue  this  area  of  oceanographic  research. 
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